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ABSTRACT 


A method is developed for finding the natural modes and frequencies 
of a vehicle with flexibly suspenc ed components. The vehicle is idealized 
as a system of one nonuniform beam and any number of flexibly suspended 
masses. A vibration analysis is first made of the beam, omitting the 
masses, by a modified Stodola method. The equations for the connected 
system are then derived using Lagrange's equations. 
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DEFINITION OF SYMBOLS 


Symbol 

A 


a ef 

B 

b ef 

C 

D 

i, 1 


i;;k) 


j,i 

k rj 


tj 


m. 

m 

m* (x) 

ra b 

M i 

n 


V q s 


Definition 
coefficient matrix 
element of coefficient matrix 
coefficient matrix 
element of coefficient matrix 
base of vehicle 
forward end of vehicle 
index of beam bending modes 
rotary inertia of beam per unit length 
moment of inertia of beam 
added moment of inertia of jth mass 
index of added mass 

rotational spring constant between main beam and jth 
attached mass 

translational spring constant between main beam and jth 
attached mass 

mass of jth component 

number of suspended masses 

mass of beam per unit length 

total mass of beam 

generalized mass of ith main beam mode 
number of single beam modes used in analysis 
general variable 
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Symbol 

T 

V 


X 

X 

me 

X 

\cj 


X 3 


X 


X. 

J 

y 


Y. (x) 
1 

y!(x) 


Y(x) 

Y' (x) 

z . 

J 

z . 

J 

0 . 

J 
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DEFINITION OF SYMBOLS (Continued) 
Definition 
kinetic energy of system 
potential energy of system 
longitudinal station coordinate 
center of mass location of beam 


X - *mc 

center of gravity location of jth component 

longitudinal station coordinate of attachment of the 
jth mass 

X. - X 
J me 

X . - X. 
mej j 

beam rigid body displacement 
displacement of beam in ith mode 
slope of displacement in ith mode 

displacement of beam in ith mode at attachment point 
of jth component 

total displacement of beam 

slope of total displacement of beam 

displacement of jth component at attachment point 

displacement of jth component at center of gravity of 
the component 

rotation of jth mass 

amplitude of ith bending mode 

main beam rotation 
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Symbol 


0) 
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DEFINITION OF SYMBOLS (Continued) 
Definition 

natural frequency of main beam in ith mode 

natural frequency of total system 

& u - 1 if i = i, 5 ii = 0 if iM 


v 
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A METHOD OF BENDING VIBRATION ANALYSIS FOR A VEHICLE 
WITH FLEXIBLY SUSPENDED COMPONENTS 

SUMMARY 


A method is developed for finding the natural modes and frequencies 
of a vehicle with flexibly suspended components. The vehicle is ideal- 
ized as a system of one nonuniform beam and any number of flexibly sus- 
pended masses. A vibration analysis is first made of the beam, omitting 
the masses, by a modified Stodola method. The equations for the con- 
nected system are then derived using Lagrange's equations. 


I. INTRODUCTION 


Space vehicles often have components of rather large mass that 
are attached to the main structure by flexible mounts. (In Apollo 
spacecraft, the lunar module is a prime example.) The natural fre- 
quency of the component may be in the frequency range of the lower 
natural frequencies of the system, in which case, the effect of these 
components on the natural mode shapes and frequencies of the vehicle 
cannot be adequately determined using a single beam model for analysis. 
Consequently, a method of analysis which allows masses to be attached 
to the main beam by flexible springs is presented. 

The vehicle is modeled as a single, continuous, nonuniform beam 
having masses attached by springs in both translation and rotation. 

An "attached mass" is assumed to be rigid. 


II. EQUATIONS OF MOTION 


The motion of the main structure is represented by a number of 
flexible body modes and two rigid body modes in one plane. (The number 
of flexible body modes is determined by the difficulty of calculation.) 
The attached components are represented as rigid bodies with flexible 
attachment, in translation and rotation, in the plane of the main body 
motion. 



The total displacement of the vehicle is given 


n 

Y(x) = ^ Y^x) + y + xcp, 
i=l 

and the slope of the displacement is denoted by 
n 

Y' ( x ) \ Yl(x) + <p. 

i=l 

The motion of the center of gravity of the attached masses is represented 
by 

z. - z. + X.6.. 

J J J J 

Lagrange’s equation 

(3l/35 a ) + - 0, s = 1,2 (n + 2 + 2n) 

results in the following equations of motion: 

At?' } + B{q 1 = (0), 

S o 

where A is the inertia matrix, B is the stiffness matrix, and 



o 



Assuming solutions of the above equations of the form 


{q s 3 = {q s 3 sin wt, 

the following equations are obtained: 


(w 2 A - B) {q g } - {0}. 


The total kinetic energy of the system is expressed by 


D n 


T = j J m' ^ Y^x) + y + X cp dx 


C i=l 

D n 


+ l J l o [I + 4>] a* 


C 1=1 


m m 


j-1 




and the total potential energy is given by 


n m 



1=1 j=l 


n 

I + * + ^ 

i=l 


z . 

3 


-2 


+ 


1 

2 


m 



k . 
r J 


n 



i=l 


+ cp) - 
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The required derivatives with time dependency eliminated follows; 


D 


D 


|r“ = T] if «* 00 Yj(x) dx + t) f X o (x) Y^(x) dx 


= f). M. 

'i 1 


— (ST/St^) = rj M ± = w 2 M i Tj i 


cjT 

Sy 


? » y m' (x) dx = y 


— (ST/S?) = m b y = -w 2 y 


D 


St 

dcp 


® cp /"* X 2 m' (x) dx = cpl^ 


(ST/Sq>) = I b cp » -w 2 I b cp 


riT 

•^r— - m. z. + m. X. 9. 
oz. 3 3 3 3 3 

3 

~r (ST/Sz.) = m.z. + m X .0 = -of m z - w 2 m X 0, 
a c j J J jj JJ jjj 
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oT 


— = m.X.z. + (m,X. +1.) e. 

de. 1 J J 3 2 23 

3 


dt ' VA + ( Vj + V 5 j 


= -co 2 m.X.z. - co 2 (m.X? + I.) 9. 
J J J 3 J 3 3 


m n 


$T— = rj. w? M. + V k .j y (tj.Y. . + y + X. 
<3r,. r l L tji /L, v ‘i lj J j 


cp) - 
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3=1 1-1 
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+ > k rjLA ^ Y ij + *» - V 

3=1 i=l 


Y*. . 
13 


z . 
3 


Y. . 
13 


m n 


V r v 

r 1 = > k . | > (r,.Y.. + y + X.cp) - z. 

cy tj ' l_j ‘i 13 J 3 W 3 


3=1 i=l 


3v 

*\ 

OCp 
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v _ r 
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r 1 

L « \L 

3=1 

i=l 

at 
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\ ' 

+ ) k . 
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'd — a 

3=1 

i=l 


Itl 
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3T = ' 

3 3=1 


k t.i 
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dV 

de. 

J 


m n 

"Z k rj [Z Vi /* 1 * 


j=l 


l=l 


Jj 


The elements of the A and B matrix are given in general form, 
is an element for each i and each j as i varies from 1 to n and j 
1 to m. 


a r - M. if i = i 
ef i 


e = i 
f = i 


a e£ "*b 


e = n -t 1 

f = n + 1 


a ef I b 


e = n + 2 
f = n + 2 


a ef = m j if J = 1 


e = n + 2 + j 
f = n + 2 + j 


= X.m. 
J J 

if j 

= j. 

= X.m. 
J J 

if j 

= i 


= (m.X. 
J J 

+ V 

if 


e = n + 2 + j 

f=n+2+m+j 

e = n + 2 + m + j 

f = n + 2 + j 

e = s n + 2+ m + j 

f=n + 2+ra + j 


The remaining elements of the A. matrix are zero. 

m m 

b f = 8 . . M, w? 4- \ k. . Y. . Y. . + ) k . Y' Y* 
ef li x l /j tj lj ij r J ij xj 


j-1 


j-1 


e 

f 


l 

i 


m 


b - = > k . Y. . 
ef Li fc J iJ 
j-1 


e 

f 


l 

n + 


There 

from 
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b 


ef 



) X. k Y. . + ) k . Y* . 
Z_j j tj xj r j ij 

j-1 j=l 


e = i 
f = n + 2 


b _ = - k . Y. . 
ef tj xj 


e = i 

f = n + 2 + j 


b 


ef 


- k . Y' . 
r 3 1*3 


e = i 

f + n 4* 2 + m ‘ 


b 


ef 


m 





Y. . 
ij 


e = n + 1 
f = i 


m 



j=l 


e = n + I 
f = n + 1 


b 


ef 


m 


j-i 


x. k 

J tj 


e * u + 1 
f = n + 2 


b 


ef 



e = n + 1 
f = n + 2 + j 


b 


ef 


= 0 


e - n + 1 
f=n+2+m+j 


ef 


m 

-V 

u 

j-1 


m 


X. k 


tj 


Y. . + 
ij 


k. . 
J-3 


y' . 


j-1 


e = n + 2 
f = i 
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b 


ef 


m 



j=l 


X. k . 
3 1 3 


e = n + 2 
f = n + 1 


b 


ef 



XT k 
J tj 



k 


e = n + 2 
f = n + 2 


b = -X. k . 
ef J tj 


e = n + 2 
f — n + 2 + j 



e = n + 2 

f = n + 2+ m-l-j 


b 


ef 


-k. . Y. . 
tj ij 


e = n + 2 + i 
f = i 



e = n + 2 + j 
f = n + 1 


b 


ef 


X. k 
3 


e = n f 2 + j 
f = n + 2 


b . = 8.. k . 
ef JJ tj 


e = n + 2 + j 
f = n + 2 + 2 


e = n + 2 + j 
f=n+2+m+j 


b 


ef 




e=n+2+m+j 
f = i 


b 


ef 


= 0 


e = n+ 2+ m+ i 
f = n + 1 
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b 


ef 



e = n+ 2+ m+ i 
f = n + 2 


b 


ef 


0 


e=n+2+m+j 
f = n + 2 + i 

Mm 


b - = 6. # k . 
ef JJ r 3 


e = n + 2+m + j 
f = n + 2 + m + j 


T, 3 solution of the equation 

(w 2 A - B) {q g } = {0) 

is acquired by transforming the equation into the standard eigenvalue 
problem. 


Ill, DISCUSSION OF RESULTS 

To evaluate the model, an analysis was made of the third stage 
(S-IVB-D, IU and payload) Saturn V dynamic test vehicle. The results 
were compared to dynamic test results recorded in references 2 and 3. 

The configuration of the vehicle is shown in figure 1. 

In the analysis, the lunar module (LEM) was idealized as a mass with 
center of gravity at the attachment point. The remainder of the vehicle 
was considered to be a heterogeneous beam. Analyses were made for engine 
ignition and engine cutoff conditions. 

The attachment spring constants used were .19017 x 10 s (newtons) 
and .1367 x 10 9 (newton meters) for translation and rotation, respectively. 
These amounts of stiffness would give LEM frequencies of 6 Hz and 10 Hz 
in translation and rotation, respectively, if the springs were attached 
to rigid structure. Four free-free bending modes of the beam were com- 
puted for input to the coupled mode analysis. Six coupled modes were 
computed, but only four are shown in this report because of lack of test 
data for comparison. 
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Dynamic test results are compared with theoretical results in fig- 
ures 2 through 11. The relative displacements obtained from the dynamic 
tests are shown for three different locations on the LEM. Only the rela- 
tive displacement at the attachment location is shown for the theoretical 
analysis. The relative displacements are normalized to unity at the place 
of maximum amplitude of the beam. (The relative displacement of the LEM 
may be larger -han unity.) The corresponding natural frequencies are 
given on each figure. 

The correlation between dynamic test results and theoretical results 
is considered acceptable. A better correlation could be expected if the 
mass characteristics of the dynamic test vehicle and the spring constants 
of the LEM were better defined. 


CONCLUSIONS 


This method of analysis is particularly useful when the natural fre- 
quency of some component is in the frequency range of the lower natural 
frequencies of the system. Both the natural frequencies and the mode 
shapes are more adequately determined using this method of analysis. 

The accuracy of results obtained from this method is highly dependent 
on a knowledge of the effective spring constants of the corrections between 
components and vehicle structure. These effective spring constants are 
often difficult to calculate by analytical means. Hence, early component 
testing to determine these spring rates would prove extremely advantageous 
when this type of analysis is needed. 
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FIG. 3. DEFLECTS CURVE 
(IGNITION FUEL CONDITION) 
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FIG. 4. DEFLECTION CURVE 
(IGNITION FUEL CONDITION) 
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FIG. 5. DEFLECTION CURVE 
(IGNITION FUEL CONDITION) 
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FIG. 6. DEFLECTION CURVE 
(IGNITION FUEL CONDITION) 
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FIG. 7. DEFLECTION CURVE 
(IGNITION FUEL CONDITION) 
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FIG. 8. DEFLECTION CURVE 
(CUTOFF FUEL CONDITION) 
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FIG. 9. DEFLECTION CURVE 
(CUTOFF FUEL CONDITION) 
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FIG. 10. DEFLECTION CURVE 
(CUTOFF FUEL CONDITION) 
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FIG. 11 DEFLECTION CURVE 
(CUTOFF FUEL CONDITION) 
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